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The distribution of gluon fields in hadrons is of fundamental interest in QCD. Using lattice 
QCD we have observed the formation of gluon flux tubes within three quark (baryon) and quark 
plus antiquark (meson) systems for a wide variety of spatial distributions of the color sources. In 
particular we have investigated three quark configurations where two of the quarks are close together 
and the third quark is some distance away, which approximates a quark plus diquark string. We 
find that the string tension of the quark-diquark string is the same as that of the quark-antiquark 
string on the same lattice. We also compare the longitudinal and transverse profiles of the gluon 
flux tubes for both sets of strings, and find them to be of similar radii and to have similar vacuum 
suppression. 
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I. INTRODUCTION 

There has recently been a renewal of interest in the 
properties of diquarks in hadronic systems, as they may 
play an important role in the existence of exotic states, 
such as the + , or in explaining the scarcity of such 
exotics 1 . While string-type models consisting of a 
quark plus a diquark have been studied using various 
analytic techniques (see [2] for a review), there have 
only been a few studies of this type of system on the 
lattice. These studies have investigated the mass of di- 
quarks [31 HI |3J O |7] , and, more recently, the nature of 
diquark correlations [8] . Recently we investigated the for- 
mation of flux tubes in static baryon systems on the lat- 
tice using a high statistics approach which enabled us to 
observe correlations between the vacuum action density 
and the quark positions in a gauge independent manner 
[llj In that work the three quarks were positioned ap- 
proximately equidistant, and a Y-shaped flux-tube was 
observed to form at large inter-quark distances. In this 
work we extend our study of three quark systems to the 
case where two of the quarks are close together and the 
third is some distance away. 

In QCD, two quarks close together, a diquark, can 
transform either according to the conjugate representa- 
tion (3) or the sextet (6) representation of SU(3). The 
color hyperfine interaction then leads to attraction in 
the spin singlet, scalar diquark channel, while the spin 
triplet, axial vector diquark is disfavoured. Hence low- 
lying diquarks should have positive parity and belong to 
the color 3 representation, and so have many properties 
similar to an antiquark. In lattice QCD this should lead 
to the formation of quark-diquark flux tubes with sim- 
ilar physical characteristics to those of quark-antiquark 
flux tubes. In particular we would expect the long range 



linear part of the quark-diquark potential to have the 
same slope as that of the quark-antiquark potential, cor- 
responding to the flux tubes having the same energy den- 
sity, and we would expect the flux tubes to have similar 
transverse size. In this work we investigate whether these 
similarities do indeed hold. 



II. FLUX TUBES ON THE LATTICE 

In order to study flux-tubes on the lattice, we begin 
with the standard approach of connecting static quark 
(and antiquark) propagators with spatial- link paths in a 
gauge invariant manner. We use APE-smeared spatial- 
link paths to propagate the quarks from a common origin 
to their spatial positions as illustrated in Fig. [T] In earlier 
work we saw that after approximately 30 APE smearing 
steps we had obtained optimal overlap with the ground 
state. The static quark propagators are constructed from 
time directed link products at fixed spatial coordinate, 
Yii Ut{x, U), using the untouched 'thin' links of the gauge 
configuration. In principle, the ground state is isolated 
after sufficient time evolution. Finally smeared- link spa- 
tial paths propagate the quarks back to the common spa- 
tial origin. 

The three-quark Wilson loop is thus defined as: 

W 3Q = }-e abc 6 a ' b ' c ' Uf Uf' Uf , (1) 
where Uj is a staple made of path-ordered link variables 

U 3 = Fexp jf dXpA»(x) \ , (2) 

and Tj is the path along a given staple as shown in Fig. [T] 
In contrast the quark-antiquark Wilson loop is given by 
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FIG. 1: Gauge-link paths or "staples," Ui, U2 and U3, forming 
a three-quark Wilson loop with the quarks located at ri, r*2 
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denote colour anti-symmetrisation at 



the source and sink respectively, while r indicates evolution 
of the three-quark system in Euclidean time. 
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In this work we have used 300 quenched QCD gauge 
field configurations created with the 0(a 2 )-mean-field 
improved Luscher-Weisz plaquette plus rectangle gauge 
action [IB] on 16 3 x 32 lattices with the long dimension 
being the x direction, making the spatial volume 16 2 x 32. 
Two hundred of these configurations were at at /3 = 4.60 
(as in our previous work) and one hundred at f3 = 4.80, 
to investigate the use of a finer lattice. These values of j3 
give lattice spacings a of 0.123 fm and 0.0945 fm respec- 
tively. 

We use lattice symmetries to improve the signal to 
noise ratio of our measurements. These include trans- 
lational invariance (any point on the lattice can be taken 
as the origin), reflection in the x plane and 90° rotational 
symmetry about the x-axis. The advantage of this ap- 
proach is that we do not have to perform any gauge fixing 
to find a signal in our flux distributions. 

We characterise the gluon-field fluctuations in our 
configurations using the gauge-invariant action density 
S(y, t) observed at spatial coordinate y and Euclidean 
time t relative to the origin of the Wilson loop. We calcu- 
late the action density using the highly-improved C(a 4 ) 
three-loop improved lattice field-strength tensor [T^] on 
four-sweep APE-smeared gauge links. While the use of 
this highly-improved action suppresses correlations close 
to the quark positions, it gives good resolution of the 
flux-tube correlations we are interested in. 

Defining the quark positions as r*, relative to the ori- 
gin of the Wilson loop, and denoting the Euclidean time 
extent of the loop by r, we evaluate the following corre- 
lation functions 



CzQ{y;ri,r 2 ,r3;T) 



(^3Q(ri,r 2 ,r 3 ;r)5(y,r/2)) 

(W 3Q (r))(S(y,T/2)) 
{W QQ {fi,r 2 ;r) S(y,r/2)) 
(W QQ (r))(S(y,T/2)) ' 



(4) 
(5) 
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FIG. 2: Projection of the T-shape path on the x-y plane. 

the product of two staples 

W QQ = 5 ab 5 a ' b ' Uf (uf\. (3) 

The three-quark configurations we use to approximate 
a quark-diquark string are T-shapes, with the origin at 
the junction of the T. Two quarks are positioned one 
lattice step in opposite directions from the origin (ap- 
proximating the diquark) , and the third is placed from 1 
to 12 lattice steps in an orthogonal direction, as shown 
in Fig. [2] 



where (• • • ) denotes averaging over configurations and 
lattice symmetries. These correlate the quark positions, 
via the Wilson loops, with the gauge-field action in a 
gauge invariant manner. For fixed quark positions and 
Euclidean time, the correlation functions are scalar fields 
in three dimensions. For values of y well away from the 
quark positions fi, there are no correlations and C — ► 1. 
Also the correlators are positive definite, eliminating any 
sign ambiguity on whether vacuum field fluctuations are 
enhanced or suppressed in the presence of static quarks. 
We find that C is generally less than 1, signaling the 
expulsion of vacuum fluctuations from the interior of 
hadrons. 

In Figs. [3] and [3] we show examples of the expulsion of 
vacuum fluctuations and the formation of flux-tubes for 
our quark-diquark and quark-antiquark configurations. 



FIG. 3: Expulsion of gluon-field fluctuations from the region 
of static quark sources illustrated by the spheres. An iso- 
surface of C(y) is illustrated by the translucent surface. A 
surface plot (or rubber sheet) describes the values of C(y) for 
y in the quark plane, (j/i, 1/2, 0). 





FIG. 5: Effective static quark potential for 30-sweep smeared 
quark-diquark sources for /3 = 4.6 (top) and (3 = 4.8 (bot- 
tom). From bottom up, the lines correspond to quark-diquark 
separation increasing from 1 to 12 lattice spacings. 



FIG. 4: Formation of quark-antiquark flux-tube. Details are 
as in [3] 

III. EFFECTIVE POTENTIALS 

In this section we extract the long range portion of 
the potential for our quark-antiquark and quark-diquark 
flux-tubes. According to QCD, the effective potentials 
should have the same slope, so long as the APE smear- 
ing in the spatial directions has smoothed the gluon fluc- 
tuations sufficiently to isolate the ground state and the 
propagation of the Wilson loops in the time direction is 
long enough for any excited states to decay. The effec- 
tive potential is obtained from the Wilson loops in the 
standard manner: 

■ v « T »- b '(^Ti))- (6) 

As shown in Fig. [5] and [6] we obtain stable plateaus for 
the potentials as a function of r. The statistical uncer- 



tainties are estimated using the jackknife method [14] . 

The quark-antiquark potential has the well-known 
form 

a 

Vqq{t) =V -- + a QQ r, (7) 

where a is the string tension. The three quark potential 
is [HUTU] 

(8) 

j<k J 

where Cf = 4/3 and L(r) is a length linking the quarks. 
As shown in our earlier work [11] , L(r) is given by the 
minimum length of string that connects the three quarks, 
or the sum of distances from the quarks to the Fermat 
(or Steiner) point. QCD suggests that the two string 
tensions ctqq and a^Q are equal. In Fig. |7] we plot 
the extracted effective potentials for the quark-diquark 
and quark-antiquark flux-tubes at each of the values 
of (3 for our gauge configurations. The plots in fermi 
show that the QCD prediction is confirmed at both val- 
ues of p. Converting length measurements from lattice 




FIG. 6: Effective static quark potential for 30-sweep smeared 
quark-antiquark sources for f3 = 4.6 (top) and j3 = 4.8 (bot- 
tom). From bottom up, the lines correspond to quark-diquark 
separation increasing from 1 to 12 lattice spacings. 



units to fermi we obtain the quark-diquark string tension 
03 q = 0.97±0.01 GeV fm _1 , which is in excellent agree- 
ment with the quark-antiquark string tension <Jqq = 0.98 
GeV fm _1 , defining the lattice spacing. 



IV. FLUX-TUBE PROFILES 

We can gain further insight into the properties of 
the flux tubes by examining their profiles close to the 
quark. We study the values of the correlators CzQ{y) 
and CgQ^y) where y = (yi,y2,0) is constrained to the 
plane of the color sources, and the origin is at the po- 
sition of either the antiquark or the join of the T. The 
quark is then at the position (£, 0, 0) where £ varies from 
1 to 12 lattice steps. 

First we examine the longitudinal profiles of both 
quark-diquark and quark-antiquark flux-tubes along the 
line (y) — (x, 0, 0) in Fig. [8] As expected, the vacuum 
expulsion close to the diquark is stronger than in the 
vicinity of the antiquark. However, near the quark the 
two flux tubes show very similar profiles. Similar results 




FIG. 7: Comparison of quark-antiquark and quark-diquark 
effective potentials for /3 = 4.6 (top) and j3 = 4.8 (bottom). 
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FIG. 8: Comparison of longitudinal flux tube profiles for f3 
4.8 at longitudinal separation of 12 lattice units. 



are seen at (3 = 4.6. 

Next we examine the transverse profiles along a line 
orthogonal to the midpoint of the flux tube, it. along 
(C/2, y, 0) for £ even, or along ((£ + l)/2, y, 0) for £ odd. 
In Fig. [9] we show profiles of both quark-diquark and 
quark-antiquark flux-tubes for £ = 12. We find that as 
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FIG. 9: Transverse profiles for quark-diquark (red lines) and 
quark-antiquark (black lines) flux tubes at f) = 4.6 (solid 
lines) and /3 — 4.8 (dashed lines). 



TABLE I: Values of the fit parameters of the function I — 
A exp(— y 2 /r 2 ) to the transverse profiles of the flux tubes, r 
is reported in lattice units (LU) and fm. The last column 
indicates the area under C(y) = 1 for the fitted curve in units 
of relative-depth times LU. 



Flux Tube 


P 


A 


r (LU) 


r (fm) 


Ar^/n 


QQ 


4.6 


0.0511(f) 


2.89(3) 


0.355(3) 


0.261(3) 


QQQ 


4.6 


0.0510(f) 


2.91(2) 


0.358(3) 


0.263(3) 


QQ 


4.8 


0.0243(3) 


3.17(4) 


0.299(4) 


0.136(3) 


QQQ 


4.8 


0.0241(3) 


3.21(3) 


0.303(3) 


0.137(3) 



V. CONCLUSIONS 

Wc have directly compared gluon flux-tubes for quark 
plus antiquark and three quark systems. In the three 
quark systems we kept two quarks close together (two 
lattice units separation), so that the system would ap- 
proximate a quark-diquark string. We found that the 
string tension in the quark-diquark string was the same 
as for the quark-antiquark string. In addition wc com- 
pared the vacuum expulsion in both sets of flux-tubes. 
We found that, in the vicinity of the quark, there was 
no measurable difference between the transverse profiles 
of the quark-diquark flux-tubes and the quark-antiquark 
flux-tubes. Also the longitudinal profiles of both sets of 
flux-tubes were very similar. 

These findings confirm the expectation from QCD that 
a diquark has many properties in common with an an- 
tiquark. In particular the long range color interaction 
between a diquark and a quark is seen to be the same as 
that between an antiquark and a quark. This result is 
interesting in that it is obtained in the quenched approx- 
imation, where the color hyperfine interaction should be 
small. This implies that the APE smearing and prop- 
agation in Euclidean time we have performed has been 
sufficient for decuplet baryon states to decay. It would 
be interesting to repeat this work with dynamical quarks, 
where variation in the strength of the color hyperfine in- 
teraction could be investigated. This is potentially of 
great importance to phenomenological models of hadron 
structure. 



long as £ is larger than one third of the total length of 
the quark-antiquark system, the transverse profiles are 
close to identical. 

The transverse profile of the flux-tube is fitted well by 
a Gaussian function C(y) = 1 — Ae~ v l r . The fit enables 
us to estimate the radius (r) and area (Ar^/n) of the flux- 
tube as well as its depth. The fit parameters for our flux 
tubes are given in Table |T] Again, the fits show that for 
long enough flux tubes the transverse profiles of quark- 
diquark and quark-antiquark flux-tubes are statistically 
identical. 
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